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We present continuum extrapolated lattice QCD results for up to fourth-order diagonal and off- 
diagonal quark number susceptibilities in the high temperature region of 300 — 700 MeV. Lattice 
QCD calculations are performed using 2-1-1 flavors of highly improved staggered quarks with nearly 
physical quark masses and at four different lattice spacings. Comparisons of our results with recent 
weak coupling calculations yield reasonably good agreements for the entire temperature range. 


I. INTRODUCTION 

At high temperatures strongly interacting matter un¬ 
dergoes a deconfinement transition, where thermody¬ 
namics can be described in terms of quarks and gluons 
[in. Lattice QCD studies [ll-[3as well as results from 
heavy-ion collision experiments U suggest that at least 
up to temperatures couple of times larger than the transi¬ 
tion temperature quark-gluon plasma (QGP) may behave 
as a strongly coupled liquid. The asymptotic freedom of 
QCD guarantees that at sufficiently high temperatures 
QGP becomes weakly coupled and should be described 
by weak coupling expansion results. However, nonpertur- 
bative effects could still remain important even at very 
high temperatures due to the infrared problems arising 
from the chromomagnetic sector [^. Thus, quantitative 
validations of weak coupling QGD calculations against 
fully nonperturbative lattice QGD results are necessary 
to ascertain the temperature range where the strongly 
coupled QGP liquid goes over to a weakly coupled quark- 
gluon gas. 

Fluctuations of and correlations among the conserved 
charges are known to be sensitive probes of deconfine¬ 
ment and are also suitable for testing the weakly or 
strongly coupled nature of QGP. The study of fluctua¬ 
tions and correlations of conserved charges on the lattice 
was initiated some time ago At low temperature 

fluctuations and correlations of conserved charges can be 
understood in terms of an uncorrelated hadron gas [ 13 - 
[i3. Deconfinement is signaled by a sudden breakdown 
of such a hadronic description in the vicinity of the QGD 
transition temperature [12, [la| ■ 

Fluctuations of and correlations among the conserved 
charges are defined through the appropriate derivatives 
of the pressure, p with respect to the chemical potentials 
associated with the corresponding conserved charges. 
In 2-1-1 flavor QCD there are three chemical potentials 
corresponding to baryon number, electric charge and 
strangeness. Since at high temperatures the relevant de¬ 
grees of freedom are quarks and gluons, it is natural to 
use the flavor chemical potentials corresponding to the 
u (up), d (down) and s (strange) quark numbers instead 
of the three conserved charge chemical potentials. In the 


flavor basis the fluctuations and correlations of charges 
get mapped into the diagonal (xl) and the off-diagonal 
(Xnm) quark number susceptibilities 

5” (p/T^) , (p/T^) 

d{pqlTf' d{pqlTf{pq,lTT' ^ ’ 

Here, pq and pqi are the chemical potentials correspond¬ 
ing to quark flavor q and g', with g, g' = u, d, s, and n, m 
denote the number of derivatives taken with respect to 
the quark flavors. 

At sufficiently high temperatures the diagonal and 
off-diagonal quark number susceptibilities should be de¬ 
scribed by weak coupling expansion results. Second- 
order quark number susceptibilities have been studied 
in weak coupling expansion for some time [isl - l^ . Re¬ 
cently also weak coupling results for the fourth-order 
fluctuations and correlations have been presented [23 - 
[27j| . The validity of these weak coupling expansion re¬ 
sults for the second-order diagonal susceptibilities X 2 has 
been tested against the accurate continuum extrapolated 
lattice QGD results in Ref. [l^]. Here we extend this 
previous study to perform accurate continuum extrapo¬ 
lated lattice QGD calculations of fourth-order diagonal 
quark number susceptibility x\ and subsequently com¬ 
pare them with weak coupling perturbative QGD results 
to validate the weak coupling regime of QGP. Diagram- 
matically, the off-diagonal quark number susceptibilities 
only consist of quark-line disconnected diagrams and cor¬ 
rections to the tree level generically start at higher orders 
of the coupling. Thus, these off-diagonal susceptibilities 
provide more stringent tests of the weak coupling regime 
of QGP. Our present study also includes lattice calcula¬ 
tions and their comparisons with weak coupling results 
of up to fourth-order off-diagonal quark number suscep¬ 
tibilities. 

The rest of the paper is organized as follows. In Sec. H 
we describe the details of our lattice QGD calculations. 
Sec. HI is dedicated to the discussion of fourth order 
diagonal quark number susceptibilities, including the de¬ 
tails of continuum extrapolations and comparison with 
weak coupling approaches. In Sec. IV we present our re¬ 
sults for the off-diagonal quark number susceptibilities. 
Finally, Sec. V contains our conclusions. 
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II. DETAILS OF THE LATTICE SIMULATIONS 


We performed lattice calculations in 2+1 flavor QCD 
using the highly improved staggered quark (HISQ) ac¬ 
tion [ 2 ^. Lattice sizes were chosen to be x Nt 
with Nt = 6 , 8 , 10 and 12 and a hxed aspect ratio of 
Na/Nr = 4. The gauge configurations used in this study 
were generated by the HotQCD Collaboration using the 
physical value of the strange quark mass rrig and degener¬ 
ate up and down quark masses = ma = nii = mg /20 
[ 33 , [3l|. The latter corresponds to a pion mass of 
161 MeV in the continuum limit [s^. For each value 
of Nt the temperature was varied by varying the lat¬ 
tice spacing a or equivalently the lattice gauge coupling 
P = 10/g^. The lattice spacing has been hxed using the 
ri scale dehned in terms of the static quark potential 


r 


2 


dr 


= 1 . 

r=ri 


( 2 ) 


We used the parametrization of a/ri from Ref. (sTj. As 
in Ref. [s^ we use ri = 0.3106 fm to convert to phys¬ 
ical units. To extend the temperature coverage in our 
calculations additional gauge conhgurations have been 
generated for /3 = 10/g^ = 8.0, 8.2 and 8.4 on 32^ x 8, 
40^ X 10 and 48^ x 12 lattices. As in Refs. we use 

the Rational Hybrid Monte Carlo algorithm [^. About 
2000 — 4000 gauge conhgurations separated by ten molec¬ 
ular dynamic trajectories of unit length have been gener¬ 
ated for each value of the gauge coupling /3. The complete 
list of the gauge conhgurations used in this study can be 
found in the Appendix. 

The quark number susceptibilities can be expressed in 
terms of the quark matrix and its inverse, and the cor¬ 
responding formulas were given in Refs. gSm. The 
necessary operators are evaluated using the random noise 
method with unbiased estimators (see Ref. [l^ for de¬ 
tails). We used between 150 to 250 random source vec¬ 
tors to evaluate the needed operators depending on the 
lattice volume and the value of /3. We have found that 
this number of source vectors is sufficient to ensure that 
the noise due to gauge held huctuations is larger than the 
noise from the stochastic estimators. In a few cases, some 
operators have been estimated with many more random 
source vectors to ensure that there are no additional sys¬ 
tematic errors due to the limited number of source vec¬ 
tors. In the Appendix we give a detailed account for the 
number of random source vectors used in our study. The 
inversion of the quark matrix was performed in double 
precision, and the squared residual of the inversion was 
less than 10 “^®. 


III. DIAGONAL QUARK NUMBER 
SUSCEPTIBILITIES 

In this section we will discuss our results on diago¬ 
nal quark number susceptibilities. The second-order light 
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FIG. 1. The light [q = u) and strange {q = s) fourth-order 
quark number susceptibilities shown as filled and open sym¬ 
bols, respectively, for different values of the lattice temporal 
extent Nt and normalized by the results for the continuum, 
massless ideal gas limit = B/vr^. The horizontal lines 

denote the lattice ideal gas results for different values of Nt- 
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and strange quark number susceptibilities have been dis¬ 
cussed in Ref. [ 2 ^ in great detail using HISQ action, in¬ 
cluding the continuum extrapolation. Our results for the 
second-order quark number susceptibilities agree with 
the Hndings of Ref. [ 2 ^, and therefore we will not dis¬ 
cuss these further. In Fig. [1] we show our results for the 
light and strange fourth-order quark number susceptibil¬ 
ities X 4 (v = u, s) for Nt = 6 , 8 , 10 and 12. The cutoff 
dependence of xt is similar to that of X 2 ’*i contin¬ 
uum is approached from below as suggested by the cutoff 
effects in the lattice free theory. Similar to the case of X 2 
[ 2 ^, the size of the cutoff effects is smaller than in the 
free theory, cf. Fig. [T] The quark mass dependence is, 
however, opposite to the case of X 2 for T < 500 MeV. 
Namely, the fourth-order light quark number susceptibil¬ 
ities are below the strange quark number susceptibilities, 
while for the second-order susceptibilities xi < X 2 [Hi- 
For T > 500 MeV the quark mass effects in the fourth- 
order quark number susceptibilities are smaller than the 
statistical uncertainty. 

To perform the continuum extrapolation, we first inter¬ 
polate the lattice results as a function of the temperature 
using smoothing splines. Th e sp line interpolation is per¬ 
formed using the R-package . The errors of the spline 
are estimated using the bootstrap method. At selected 
values of the temperature in the interval of T = 340 MeV 
to T = 660 MeV, separated by 20 MeV, we performed 
continuum extrapolations using the interpolated values 
of X 4 and the corresponding bootstrap errors. Illustra¬ 
tive results for the continuum extrapolations for xl are 
shown in Fig. [21 For each temperature we performed the 
continuum extrapolation using the form motivated by the 
lattice free theory [^ . 

xl{NT)=a + b/N^ + c/Nl (3) 

We also performed extrapolations using only the data 
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FIG. 2. Continuum extrapolations of X 4 a-t three representative temperatures. In the middle panel, we show extrapolation for 
Xi and compare our results with the Nr dependence of xl obtained with p4 action [^. The solid and dashed lines correspond 
to continuum extrapolations using Eq. while the dotted lines correspond to extrapolations performed using Eq. ; see 
the main text. The filled triangle corresponds to the continuum value. All the results have been normalized by 


for Nr > 8 and setting c = 0. These extrapolations typ¬ 
ically give larger continuum values but remain compati¬ 
ble within errors with the original extrapolations. These 
fits are shown as solid and dashed lines, respectively, in 
Fig. [21 We do not find evidence for a significant l/N^ 
term in the cutoff dependence for xl- This is evident 
from the Nr dependence of our results shown in Fig. [2] 
A similar situation was observed for xt lUl- Since the 
cutoff dependence of xt qualitatively follows the lattice 
free theory expectation we also performed continuum ex¬ 
trapolation using the form 

xliNr) = e + / • - xf(4) 

where xl’is the lattice free theory result and 
X 4 ^eai = G/tt^ is the continuum free theory result for the 
fourth-order quark number susceptibility. This simple 
two-parameter fit describes the data well. The results of 
these fits are also shown in Fig. [2] as dotted lines. For the 
coefficient / we get values ranging from about 0.30 at the 
lowest temperature to about 0.43 at the highest temper¬ 
ature; i.e., cutoff effects are only 30% — 40% of those in 
the free theory. The central values of the extrapolations 
obtained from Eq. o are typically smaller than the ones 
obtained from the three-parameter fit with Eq. I®), but 
still compatible within errors. In other words, there are 
no statistically significant deviations between the differ¬ 
ent continuum extrapolations. Therefore, the extrapo¬ 
lated values obtained from three-parameter fits with Eq. 
([3]) will be used as our continuum results. Since the dif¬ 
ferent continuum extrapolations give a difference of the 
order of the statistical errors, we conservatively estimate 
the total error on our continuum results for xt us twice 
the statistical error obtained from our three-parameter 
fits. 

In Fig. 12] we also show the lattice results for xl for 
T = 500 MeV obtained with p4 action for Nr = 6,8 and 
12 [ 25 . As discussed in Ref. [l^, the continuum limit 
for p4 action is approached from abovei, and the cutoff 
dependence in this case does not follow the free theory ex¬ 
pectation; the leading-order cutoff dependence goes like 
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FIG. 3. The continuum extrapolated results for xl and xl 
normalized by the corresponding massless ideal gas value 
X4 = d/tT . 



l/Nr- As was also discussed therei, the quality of the p4 
data does not allow a reliable continuum extrapolation 
for xl- However, in the continuum limit, the results ob¬ 
tained with HISQ action and p4 action should of course 
agree. Therefore, we fit the cutoff dependence of the p4 
data using the form a + g/Nr + h/Nr and requiring that 
in the continuum limit it agrees with the HISQ result. 
This fit gives x^/d.o.f. ~ 0.1. Therefore, we can at least 
say that the p4 data are consistent with the continuum 
extrapolations performed with HISQ action. The con¬ 
tinuum extrapolated results for xl and xl are shown in 
Fig. |3|as a function of the temperature. We can see that 
xl > xl but the difference between xl and xl becomes 
smaller with increasing temperature. The continuum ex¬ 
trapolated results for xt ^re provided in Table |T| 

The deviation of xl from the ideal gas limit at the 
highest temperature studied by us is 15%. This should 
be compared to the continuum result for xl which is only 
6% below the ideal gas limit at a similar temperature [ 25 . 
It is interesting to see whether these deviations from the 
ideal gas limit can be understood in terms of weak cou- 











IV. OFF-DIAGONAL QUARK NUMBER 
SUSCEPTIBILITIES 
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FIG. 4. The continuum extrapolated result for X 4 compared 
to perturbative EQCD calculations shown as the shaded band. 
The width of the band corresponds to the variation of the 
renormalization scale from irT to dvrT. The dashed lines cor¬ 
respond to the three-loop HTL calculations evaluated for the 
renormalization scale A = AnT, 2'kT and tvT (from top to bot¬ 
tom). All results have been normalized by the corresponding 
massless ideal gas results of = h/vr^. The inset shows 

the comparison of the lattice and weak coupling calculations 
for X 2 [23 normalized by — 1- 


pling expansions. The fourth order quark number suscep¬ 
tibilities have been calculated in three-loop hard thermal 
loop (HTL) perturbation theory recently in Ref. and 
in perturbative calculations using dimensionally reduced 
electrostatic QCD (EQCD) [ 2 ^. These weak coupling 
results are compared with our continuum results for X 4 
in Fig.Sl We show the scale uncertainty of the perturba¬ 
tive results by varying the renormalization scale A from 
ttT to 47rT. The comparison of the continuum lattice re¬ 
sults for X 2 @1 with weak coupling calculations is also 
shown in Fig.|3]as an inset. The EQCD band is above the 
continuum extrapolated lattice data for X 4 , while these 
calculations give results for X 2 that agree with the lat¬ 
tice data. So there remains some tension between the 
EQCD calculations and the lattice results on second and 
fourth-order quark number susceptibilities. The three- 
loop HTL perturbative result agrees very well with our 
lattice data for the renormalization scale A = 27rT. This 
scale choice, however, overpredicts the value of X 2 - ^1" 
though the scale uncertainty in the three-loop HTL cal¬ 
culation of X 2 i® rather large, it should be noted that, 
once the renormalization scale A is fixed by comparing 
one observable with the corresponding lattice QCD re¬ 
sult all the other quantities are completely parameter-free 
predictions of HTL calculations without any uncertainty. 
Unfortunately, the lattice QCD results for xt and xt can¬ 
not be simultaneously reproduced by the three-loop HTL 
calculations with a single renormalization scale for the 
coupling. Overall the three-loop HTL perturbative re¬ 
sults agree with the lattice data within a somewhat large 
scale uncertainty. 


In this section we show our results on off-diagonal 
quark number susceptibilities and compare them with 
weak coupling calculations. The off-diagonal suscepti¬ 
bilities vanish in the infinite temperature limit, i.e., for 
the ideal quark gas. In fact, using weak coupling calcu¬ 
lations, one can show that they are related to the cou¬ 
pling of quarks to soft (static) gluons. In the language of 
EQCD, one may say that they do not receive contribu¬ 
tion from the scale 2ttT. In particular, the leading-order 
contribution to X 22 comes from the term proportional 
to TrrVj-i (the so-called plasmon term) in the expression 
of the pressure, with mu being the leading-order De¬ 
bye mass. Thus, the leading-order contribution to X 22 
is of order ^ . The leading-order contribution 

to xVi i® order ~ and is sensitive to the static 
chromoinagnetic sector. Therefore, it is nonperturbative 
and can only be calculated using lattice simulations of 
EQCD [2l|. Since the off-diagonal quark number suscep¬ 
tibilities at high temperature are mostly sensitive to soft 
gluon fields, the cutoff {Nr) dependence of these quan¬ 
tities in lattice calculations using HISQ is expected to 
be small at high temperatures. Our lattice QCD results 
for x'ti s-nd X 22 shown in Fig. [5] As expected, the 
apparent cutoff dependence of the off-diagonal quark sus¬ 
ceptibilities is considerably smaller than for the diagonal 
ones. The off-diagonal susceptibilities are small, which 
is in qualitative agreement with weak coupling expecta¬ 
tions discussed above. Since the cutoff dependence of 
the off-diagonal susceptibilities is quite small, it is possi¬ 
ble to perform a comparison with weak coupling calcula¬ 
tions using the lattice data at fixed Nr- Nevertheless, we 
performed continuum extrapolations also for off-diagonal 
quark number susceptibilities. The procedure of contin¬ 
uum extrapolation here is similar to that for diagonal 
quark number susceptibilities. We first perform an in¬ 
terpolation in the temperature and then continuum ex¬ 
trapolations at the same temperature values as before 
separated by 20 MeV. We performed continuum extrap¬ 
olations fitting the Nr dependence of the lattice data 
with a + b/Nr + c/N^- We also performed fits using only 
Nr > 8 data and setting c = 0. Both fits give result that 
agree well within the errors. In the case of X 22 '''^®'® 

also possible to set 6 = 0, i.e., perform fits to a constant 
for Nr > 8. These fits give the result with the smallest 
errors but are consistent with the above fits. In the case 
of Xii, we performed fits with c = 0 to obtain the contin¬ 
uum estimate. In addition averaging over Nr = 10 and 
Nr = 12 data for Xii gave results which are consistent 
with the above continuum estimate. We used the two- 
parameter fits and the corresponding errors for our final 
continuum estimates for X 22 and Xii- 

These continuum estimates are also shown in Fig. [S] 
The numerical values of the continuum extrapolated off- 
diagonal susceptibilities are also presented in Table ID We 
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FIG. 5. The fourth-order off-diagonal susceptibility X 22 (left) and the second order off-diagonal quark number susceptibility 
Xii (right) calculated for different values of N-r. The lattice results are compared to EQCD calculations shown as the 
shaded band corresponding to renormalization scale variation from ttT to dvrT and to three-loop HTL perturbation theory 
shown with the dashed lines for renormalization scale A = ttT, 27rT and 47rT (from top to bottom). 


also calculated Xii X 22 found that within the 
estimated errors they agree with Xii £md X 22 - 

Now let us compare our lattice results with weak cou¬ 
pling calculations in more details. In Fig. [S] we also 
show the results from three-loop HTL perturbation the¬ 
ory for X 22 using three choices of the renormalization 
scale A = ttT, 2ttT and 47rT The scale choice 

A = AttT works best for the higher temperature, while 
for a lower temperature, the choice A = 27rT seems to 
be better. Overall it is fair to say thati, within the un¬ 
certainties, the lattice and the three-loop HTL perturba¬ 
tion theory results agree. We note that the three-loop 
HTL calculation for quark number susceptibilities has 
been performed in the limit of vanishing quark masses. 
As stated above the quark mass effects are negligable for 
X 22 I uud therefore comparison of three-loop HTL results 
and the lattice results for X 22 is justified. We compare 
our results for Xii with EQCD calculations of Ref. (^ . 
In EQCD Xii can be written as 

§ (si? +'*’'=<''>)' <"> 

where y = gs and niE are the gauge coupling 

and the mass parameter of EQCD which can be cal¬ 
culated at any order in perturbation theory [s^. The 
constant /3 m is a nonperturbative constant that deter¬ 
mines the leading-order contribution to x^f. We use the 
value Pm = 0.1 determined in Ref. [^. The function 
Sxsiy) parameterizes higher order corrections and was 
calculated in Ref. [2lj. Using this as well as the next- 
to-leading-order expression for y from Ref. together 
with the next-to-next-leading-order result for gs |38j . we 
calculate the EQCD result for Xii- This result is shown 
in Fig. [S] as a band. The width of the band for Xii shown 
in Fig. [S] corresponds to the error on 6 x 3 and the scale 
variation between A = ttT and A = 47rT combined. We 
see that there is a fair agreement between the lattice re¬ 
sults and the EQCD calculation for Xii- 

We also calculated the three other fourth-order off- 


diagonal quark number susceptibilities, Xst, X 31 and x“|. 
Unfortunately, these quantities are very noisy, but within 
large errors appear to be both temperature independent 
as well as Nr independent. If we fit all the available lat¬ 
tice data in the temperature range T = 280 — 710 MeV 
with a constant, we obtain the following estimates: 

X^f = 0.00020 ± 0.00015 , (6) 

xTi = 0.00071 ± 0.00030 , (7) 

= 0.00060 ± 0.00030 . (8) 

We see that these off-diagonal susceptibilities are similar 
in magnitude to Xii which may be expected in the weak 
coupling calculations. 

V. CONCLUSIONS 

In this paper we have extended the study of lattice 
QCD calculations quark number susceptibilities at high 
temperatures using the HISQ action to include fourth- 
order diagonal as well as off-diagonal quark number sus¬ 
ceptibilities. In our calculations we use lattices with tem¬ 
poral extent Nr = 6, 8, 10 and 12, and we cover tem¬ 
peratures ranging from 300 up to 700 MeV. We have 
obtained sufficiently accurate reliable continuum extrap¬ 
olated lattice QCD results at high enough temperatures 
to perform meaningful comparisons with weak coupling 
QCD calculations. We find that the cutoff dependence of 
the fourth-order diagonal quark number susceptibilities 
x\ is similar to the cutoff dependence of the second or¬ 
der quark number susceptibility xi s-nd continuum 
extrapolations can be performed in a similar manner. 
The systematic errors of the continuum extrapolations 
are quite small, smaller than or equal to the statistical 
errors. We compared our continuum extrapolated results 
for x'i to three-loop HTL perturbation theory as well as 
to EQCD results. The result of three-loop HTL perturba¬ 
tion theory agrees very well with our continuum extrap¬ 
olated result for X 4 - However, three-loop HTL calcula- 
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tions cannot simultaneously reproduce the lattice QCD 
results for y® with the same choice of the renor¬ 

malization scale for the coupling. The EQCD result of 
Ref. [23 for Xa somewhat above our continuum lattice 
result, while it agrees with the continuum lattice results 
for X 2 [ 23 . 

In the case of the off-diagonal susceptibilities, we also 
find decent agreements between the lattice and the weak 
coupling results for Xii and X 22 - From our analysis it is 
clear that lattice results provide stringent tests for differ¬ 
ent weak coupling approaches at high temperatures. The 
agreement between the lattice and the weak coupling re¬ 
sults suggests that for T > 300 MeV quark degrees of 
freedom are weakly coupled. 

While this paper was being finalizedi. Ref. [13 ap¬ 
peared. Reference [13 performed similar lattice QCD 
calculations with a different fermionic discretization 
scheme and compared up to fourth-order diagonal and 
off-diagonal susceptibilities with the weak coupling re¬ 
sults. The lattice results and hence the conclusions of 
Ref. [13 is very similar to our present study. 


TABLE I. The continuum estimates for the diagonal and off- 
diagonal quark number susceptibilities. 


T [MeV] X 4 


340 

0.4653(48) 

360 

0.4687(41) 

380 

0.4721(34) 

400 

0.4755(30) 

420 

0.4789(29) 

440 

0.4821(28) 

460 

0.4853(29) 

480 

0.4883(29) 

500 

0.4912(30) 

520 

0.4939(30) 

540 

0.4964(29) 

560 

0.4986(28) 

580 

0.5006(27) 

600 

0.5023(25) 

620 

0.5039(24) 

640 

0.5053(24) 

660 

0.5067(26) 


Xl Xii ■ 10^ 

0.4862(58) -3.62(46) 
0.4854(43) -3.29(39) 
0.4848(40) -2.99(36) 
0.4849(44) -2.77(33) 
0.4860(47) -2.56(33) 
0.4879(47) -2.34(32) 
0.4902(47) -2.13(31) 
0.4927(45) -1.98(32) 
0.4956(44) -1.85(33) 
0.4987(43) -1.73(33) 
0.5018(43) -1.61(33) 
0.5044(41) -1.50(32) 
0.5065(38) -1.42(31) 
0.5079(34) -1.36(31) 
0.5089(30) -1.30(32) 
0.5096(28) -1.24(34) 
0.5103(32) -1.20(38) 


X2I 

0.0203(21) 

0.0188(18) 

0.0171(16) 

0.0158(15) 

0.0149(14) 

0.0141(13) 

0.0135(13) 

0.0130(14) 

0.0126(14) 

0.0123(14) 

0.0120(14) 

0.0117(13) 

0.0112(17) 

0.0107(12) 

0 . 0102 ( 12 ) 

0.0097(12) 

0.0091(13) 
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For some ensembles we used many more source vectors 
to check for systematic errors. 


TABLE 11. Parameters of the calculations for 24® x 6 lattices 


P 

T (MeV) 

primp 
^ sv.l 

^sv,s 

# configurations 

6.664 

281.1 

150 

0 

150 

0 

2990 

6.800 

320.5 

150 

0 

150 

0 

3000 

6.950 

369.8 

150 

0 

150 

0 

3000 

7.150 

445.9 

150 

0 

150 

0 

2540 

7.280 

502.5 

150 

0 

150 

0 

2720 

7.373 

546.7 

150 

0 

150 

0 

2310 

7.596 

667.2 

150 

0 

150 

0 

3000 


Appendix: Details of lattice calculations 

In this Appendix we present some details of our calcu¬ 
lations. In Tables lElIIIlllY] and El we show the gauge 
coupling P = 101 and the corresponding temperatures 
for different lattice sizes. The number of gauge configu¬ 
rations used in the analysis is shown as the last column of 
the tables. We estimate the fermion operators involving 
light (u or d) quarks using random source vectors, 
while for the operators involving strange quarks we use 
Nsv,s source vectors. We tried to improved the stochas¬ 
tic estimates by using additional source vectors for all 
operators that only require two inversions. In most cases 
these operators turn out to be the noisiest but because 
of fewer inversions are cheaper to calculate. The num¬ 
ber of additional source vectors is denoted by and 

NllpP for light and strange quark operators, respectively. 


TABLE III. Parameters of the calculations for 32® x 8 lattices 


P 

T (MeV) Ns^,i 

primp 

sv.l 

^sv,s 

Ns'v'P # configurations 

7.150 

334.4 

250 

0 

250 

0 

4020 

7.280 

376.8 

250 

0 

250 

0 

4080 

7.373 

410.0 

250 

1000 

250 

0 

3940 

7.596 

500.4 

250 

0 

250 

0 

4050 

7.825 

611.5 

250 

0 

250 

0 

3920 

8.000 

711.3 

250 

0 

250 

0 

1800 
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TABLE IV. Parameters of the calculations for 40^ x 10 lattices 


p 

T (MeV) Ns^,i 

jsrimp 

^sv,s 

771/]? # configurations 

7.373 

328.0 

300 

0 

300 

0 

4060 

7.596 

400.3 

300 

0 

300 

0 

3180 

7.825 

489.2 

300 

0 

300 

0 

3090 

8.000 

569.1 

600 

600 

200 

0 

2670 

8.200 

675.3 

200 

0 

200 

0 

1590 


TABLE V. Parameters of the calculations for 48® x 12 lattices 


P 

T (MeV) Nsv,i 

jsrimp 

^SV,S 

77//]? configurations 

7.596 

333.6 

250 

250 

250 

250 

3360 

7.825 

407.7 

250 

250 

250 

250 

3020 

8.200 

562.7 

250 

250 

250 

250 

2740 

8.400 

666.9 

250 

250 

250 

250 

2710 
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